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ON 

Abstract. Classical Bianchi-Lie, Backlund and Darboux transformations are considered. Their generalizations for the 
J_l i dynamical systems are discussed. For the transformation being the generalization of the normal shift the special class of 

dynamical systems is defined. The effective criterion for separating such systems in form of partial differential equations is 
found. 
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1. Introduction 



Backlund transformations which are well known to the specialists in integrable nonlinear equations originally first 



arose in the works of classics of differential geometry in last century. Let's remind some of their results. Consider 
2-dimensional surface S in M 3 with the standard scalar product. Let us map each point M of the surface S onto the 
, point M' of the other surface S' so that the following conditions are fulfilled 

, (1) The distance |MM'| is constant and it is equal to R; 

(2) The tangent planes r and r' of the surfaces S and S' in the points M and M' are orthogonal: Z (r, r') = 90°; 

(3) The segment MM' lies in both tangents planes r and t' of the surfaces S and S' respectively. 

Bianchi [1] showed that if S is the surface of constant negative Gaussian curvature K = —1/R 2 then in assumptions 
(l)-(3) the second surface S' also is the surface of constant negative Gaussian curvature K' = —1/R 2 . Lie [2] 
strengthened Bianchi's result having shown that construction (l)-(3) can be realized only on the surface S of constant 
negative Gaussian curvature K = —1/R 2 . 

Backlund [3] generalized the construction (l)-(3) having substituted condition (2) by more weak condition of 
O ■ constancy of the angle between the tangent planes r and r'. Such kind of transformation also is defined only on the 

surfaces of constant negative curvature and it leads to the surfaces S' of the same curvature. Darboux [4] offered 
further generalization of this construction having substituted condition (3) by condition of constancy of the angles 
between line MM' and both tangent planes r and t' . However thereby some differences there appeared: Darboux 
transformation is realized on the surfaces where some linear combination of Gaussian and mean curvatures is constant. 
The result of transformation is the surface where some other linear combination of Gaussian and mean curvatures is 
constant. 

Differential equations which are obtained in considering the transformations (l)-(3) as well as their generalizations 
became the objects of the numerous investigations. For the modern state of such investigations one can enquire 
the monographs [5] and [6]. The generalization of the above transformations for the multidimensional spaces and 
submanifolds in them is made in papers [7], [8] and [9]. The generalizations for the submanifolds in the arbitrary 
Riemannian manifold of constant sectional curvature is considered in [10] and [11]. 

The normal shift (or the Bonnet transformation) is some particular case in the general Darboux construction. In 
this case tangent planes r and t' are parallel while the segment MM' is orthogonal to them both. This is degenerate 
particular case since it can be realized on any surface imposing no limitations for its curvatures. 

In this paper we discuss the generalization of the Bianchi-Lie, Backlund and Darboux transformation for the spatially 
anisotropic case substituting the straight line segment or geodesic segment MM' by the segment of trajectory for some 
dynamical system. Let's consider the following dynamical system in W 1 

(1.1) r = F(r,r) 
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This is the natural second order dynamical system defined by the force F. For each point M on some submanifold 
S C R n let us consider the particle starting from it with some initial velocity v. In the end of some time interval 
(same for all particles) these particles form some other submanifold S' C M. n . In most general form the problem may 
be stated as follows: what kind of limitations for the dynamical system (1.1) itself and for the choice of submanifold 
S and initial velocities of particles on it arise if one require the angles defining the mutual arrangement of tangent 
spaces r, r' and particle trajectories to be constant. Such problem has a lot of possible specializations one of which 
leading to the meaningful results is considered below. 

2. The normal shift along the dynamical system. 

Let's consider the dynamical system (1.1) in Euclidean space W 1 with the standard scalar product. Let S be the 
submanifold of codimcnsion 1 in M™ and let n be the vector field of unit normal vectors on S. We direct the initial 
velocity of particles along n defining the modulus of velocity as some smooth function v = v(M) on the submanifold 
S. Then the dynamical system (1.1) produces the family of submanifolds S t 

(2.1) ft : S — S t 

together with the diffeomorphisms ft binding S with St- 

Definition 1. Each transformation / = f t of the family (2.1) is called the normal shift along the dynamical system 
(1.1) if each trajectory of (1.1) crosses each submanifold St along its normal vector n. 

Definition 2. Dynamical system (1.1) is called the dynamical system accepting the normal shift of submanifolds 
of codimcnsion 1 if for any submanifold S of codimcnsion 1 there is the function v = v(M) on S such that the 
transformation (2.1) defined by the system (1.1) and the initial velocity function |v(M)| = v(M) is the transformation 
of normal shift. 

Note that the transformation (2.1) may implement the normal shift for some particular submanifolds even when 
they do not satisfy the definition 2. Dynamical systems accepting the normal shift for arbitrary submanifold form the 
special class of dynamical systems narrow enough to be described in much details. In the following two sections we 
consider such dynamical systems in R 2 and derive the partial differential equation for the force function F of them. 

3. Dynamical systems in R 2 accepting the normal shift. 

Let's consider the second order dynamical system (1.1) in Euclidean space M 2 with the standard scalar product. 
Phase space of the system (1.1) is four-dimensional in this case. The locus of points where v = r = is the two- 
dimensional plane in it. Everywhere out of this locus we define the unit vector 

(3.1) N = N(v) = ^ 

and the unit vector M(v) perpendicular to N(v). Using (3.1) the right hand side of (1.1) can be rewritten as follows 

(3.2) F(r, v) = A(r, v)N(v) + B(r, v)M(v) 
Cartesian components of N(v) and M(v) satisfy the following differential equations 

dN k MiM k dM k M,N k 



(3.3) 



dv l Ivl dv l 



Gradients of the functions A(r,v) and B(r,v) in (3.2) also can be expressed in terms of components of N(v) and 
M(v) 

dA dA 

(3.4) — t = ai N t + a 2 Mi — t = a 3 N, + a A M t 

(3.5) — t = p 1 N l + p 2 M, — t = foNi + /3 4 M 4 
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For the dynamical system (1.1) with the right hand side (3.2) we consider the Cauchy problem with the initial data 
depending on some scalar parameter s 

(3.6) r(t,«)| t=0 = r(s) 8r(t, s)\ t=0 = v(s) 

Let t(£, s) be the derivative of r(t, s) by the parameter s. Differentiating (1.1) we derive the following equation for 
the components of the vector r(i, s) 

Using the expression (3.2), the differential equations (3.3) and the formulae (3.4) and (3.5) we my write the partial 
derivatives in (3.7) in the following form 

(3.8) — = (aiJVi + a 2 M 4 ) N k + (ftiVj + foMi) M k 



- = a 3 NiN k + (04 - -^j MiN k + 



+ foNiM k +((3 4 + —) MiM 



(3.9) dV 



Time derivative for the vector v is defined by the equation (1.1). For the length of this vector then we have 
(3.10) d t \v\=A a t (|v|- 1 )=-A 

For the derivatives of A and B in (3.4) and (3.5) we obtain 

dtA = ai|v + a 3 A + a 4 B 
(3 ' H) d t £ = /3i|v|+/3 3 A + /? 4 B 

Taking into account (3.4) and (3.5) from (3.3) and (3.10) we derive the formulae for the time derivatives of N and M 

(3.12) a ,N = ^i 6,M = -^ 

M |v| 

Time dynamics of the vector r due to (3.7) and the relationships (3.12) determine the dynamics of the scalar 
products (r, N) and (r, M). Let's denote them as follows 

(3.13) (t,N)= V {tM)=4> 
Differentiating (3.13) and taking into account (3.12) we obtain 

(3.14) l l l 

(d t T,M) =d^ + — <p 

Differentiating the left hand sides in (3.14) we get 

d t (dfT, N) = (d tt T,N) + (d t r,d t N) = 

(3.15) B B 2 

| V | |v| 2 

dt (d t r, M) = (d tt r, M) + (d t r, 8 t M) = 

(3.16) b B 2 

= (^t,M)--%,+ --^/> 
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Differentiating the right hand sides in the same equations (3.14) and taking into account (3.10), (3.11) and (3.12) we 
have 

d t (d t r, N) = dttip - rid t i> + r4>ip~ 
(3-17) 1 1 11 

d t (d t r, M) = d tt ip + %> - — p+ 

(3.18) A B 

+ pup + p 3 ip + pi <p 

The second derivative of the vector r in (3.15) and (3.16) can be calculated on the base of (3.7), (3.8) and (3.9). Then 
we get 

(d U T, N) = anp + a 2 ip + a 3 (d t <p - p r V ) + 

(3.19) , x , V |V| J 

( B \ / B 

+ I a 4 — 1 — T ) ' ( 9 t tp + -r—rif 



(d tt T, M) = p x ip + p 2 ip + fo [ dtif - T^rV ) + 



(3 - 20) A\ / , B 



|v| 

As a result of equating (3.15) with (3.17) and (3.16) with (3.18) after substituting (3.19) and (3.20) for <p and ip we 
obtain 



B BA Aip Bib 

dtt^P ~ -r-rdtip + |-|2 V - Pli> ~ 03-r-r ~ ^tt = 
|v| |v| J |v| |v| 

= amp + a 2 tp + a 3 (dtp - 7-7 V'] + a 4 (dtip + ^<p 



(3.21) 

+ A d t <p - + Pup + P 3 ^j + Pa^j 
|v| |v| 2 |v| |v| 

(3.22) = Pi<P + M + (fo- ■ (dttp- + 

+ a 4 (p 4 + |^ • (dftp + -^ip 

Let the dynamical system (1.1) with the force field (3.2) be accepting the normal shift in R 2 . Submanifolds of 
codimension 1 in M 2 are the plane curves. It is convenient to use the natural parameter on them measuring the arc 
length referenced to some fixed point on the curve. Let the parameter s in the Cauchy problem initial data (3.6) do 
coincide with the natural parameter on the curve. Then for the transformations of the normal shift these initial data 
may be rewritten as follows 

(3.23) r(t,«)| t=0 =r(«) d t r(t, s)| t=0 = v(s)n(s) 

The derivative d s r(s) = r(s) is the unit tangent vector for the curve while n(s) is the unit normal vector for it. 

From (3.23), (3.14) and from the orthogonality of the vectors r(s) and n(s) we derive the following initial data for 
the function <p introduced in (3.13) 

(3.24) ip\ t=0 = d t ip\ t=0 = d s v(s) + ^ 
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By the proper choice of the function v(s) in (3.23) and (3.24) (see also the definition 2) one can vanish the right hand 
side of the second initial data statement in (3.24). Such proper choice is defined by the following differential equation 

The right hand side of (3.25) as well as the function v(s) itself then depend on the form of the curve. Once the choice 

(3.25) for v(s) is made the initial data (3.14) take the form 

(3.26) <p\ t=0 = d t <p\ t=0 = 

Now let's come back to the differential equations (3.21) and (3.22). They are linear with respect to ip and tp and 
altogether they form the complete system of the differential equations. As for the initial data (3.26) they are not 
enough to determine the Cauchy problem for this system. However one can make the special choice of the functions 
A and B in (3.2) such that the coefficients of tp and dtip in (3-21) vanish making (3.21) the separate equation of the 
second order with respect to the function ip. Because of the linearity the Cauchy problem (3.26) for this equation has 
the unique solution being identically zero. In this case vectors r and N are perpendicular for all s and t which is just 
the condition of normal shift. The proper choice of A and B leading to this case is described by the following theorem. 

Theorem 1. Dynamical system (1.1) with the right hand side of the form (3.2) is the system accepting the normal 
shift in R 2 if and only if the following relationships 

(3.27) B = -|v|a 4 

BA A B B 

(3.28) /? 4 =a2 _ a3 

|v| 2 |v| |v| |v| 

are fulfilled. Here ct\, a.i, ctz, oti and [3i, (32, (3%, (3a are the coefficients in (3.4) and (3.5). 

4. Some examples in R 2 . 

Because of the derivatives in (3.4) and (3.5) the equations (3.27) and (3.28) are the system of the partial differential 
equations with respect to A and B defining in turn the force field F according to (3.2). In order to make these 
equations more explicit one should use the explicit form of the vectors N and M 



v 



M= FT 
v 



-v 2 



For the function B then from (3.27) we may obtain its expression via the function A 
t A^ d 2dA l dA 

Further substitution of (4.1) into (3.5) and (3.28) let us bring the system of equation (3.27) and (3.28) to the form 
of one equation with respect to A. We do not write this equation here because it is very huge. But in place of it we 
consider some examples when this equation can be simplified to rather observable size. 

Case 1. Spatially homogeneous force field directed along the vector v of velocity. Function A in this case does 
not depend on r and B is equal to zero identically. From (4.1) then we obtain that A depends only on the modulus 
of velocity A = A(\v\). Therefore ct2 — and the equation (3.28) become the identity. Note that this case is trivial 
from geometrical point of view since trajectories of particles are straight lines and associated transformation in (2.1) 
coincides with the classical normal shift. 

Case 2. Spatially homogeneous but anisotropic force field. Both functions A and B do not depend on r. Let's 
denote v = v the modulus of velocity and denote via 9 the angle between v and some fixed direction in space. Then 
B = —dgA = —Ag and the equation (3.28) is written as follows 

(4.2) AAg - vAA 9v + AgAgg = -vA e A v 
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Equation (4.2) has the particular solution with the separated variables A = A(v) cos(6>). Force field then is of the form 
(4.3) F = A(v)N cos(8) + A(v)M sin(0) 

The modulus of force here |F| = A(v) is some arbitrary function of v. The direction of force F form with the fixed 
direction in space the angle 29 twice as greater than the angle between v and that direction (see fig. 1). 



Place for figA. Place for fig-2. 



Case 3. Spacially unhomogeneous force field with the central point. Here it is convenient to introduce new variables 
p, 7, v and 9 according to the fig. 2 



( 4 - 4 ) 



r 1 = pcos(7) v 1 — wcos(7 + 9) 

r 2 = /9sin(7) v 2 = wsin(7 + 9) 



From (4.4) it is not difficult to derive the following formulae for different functions in the equations (3.27) and (3.28) 

a 2 = - sm(9)A p + C ° & ^ (A 1 - A e ) a 4 = -A e a 3 = A v 
p v 

p 1 ^cos(9)B p + ^l(B J -B e ) 0A = ^B e fc = B v 

Let's substitute these expressions into (3.27) and (3.28) and impose one additional condition A 7 = (the condition of 
isotropy for all rays coming out from the central point). As a result we obtain the following equations 



B = -A. 



AA e ... sin(6>) AA 6v A e A ee 
(A a - + co ^)Ae p A ee H ~— = 

■ //,\ a cos(9) AgA v 

= - sm(9)A p K -^A e + — 2 — H 

p v 

The solution of the equations (4.5) with separated variables here has the following form 

A = ^lco S (9) B=^± S in(9) 
P P 

where A(v) is some arbitrary function of modulus of velocity v. Corresponding dynamical system has the following 
force field 

p = AHNco S (fl)+^)M S in(fl) 
P 

(see fig. 2). Note, that both dynamical systems with force fields (4.3) and (4.6) are integrable via quadratures. 
Moreover for some special choice of A(v) they are explicitly integrable. For instance when A(v) = const the trajectories 
of the system (4.3) are the cycloids. 
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